Estimating reaction rates and uncertainties for primordial nucleosynthesis 
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We present a Monte Carlo method for direct incorporation of nuclear inputs in primordial 
nucleosynthesis calculations. This method is intended to remedy shortcomings of current error 
estimation, by eliminating intermediate data evaluations and working directly with experimental 
data, allowing error estimation based solely on published experimental uncertainties. This technique 
also allows simple incorporation of new data and reduction of errors with the introduction of more 
precise data. Application of our method indicates that previous error estimates on the calculated 
abundances were too large by as much as a factor of three. Since uncertainties in the BBN calculation 
currently dominate inferences drawn from light-element abundances, the re-estimated errors have 
important consequences for cosmic baryon density, neutrino physics, and lithium depletion in halo 
stars. Our direct method allows detailed discussion of the status of the nuclear inputs, by identifying 
clearly the places where improved cross section measurements would be most useful. 
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I. INTRODUCTION 

Big-bang nucleosynthesis (BBN) is an important com- 
ponent of the hot big-bang cosmology. It provides a 
direct probe of events less than one second after the 
big bang, as well as key evidence for the existence of 
non-baryonic dark matter. The success of big-bang 
nucleosynthesis theory is indicated by the narrow range 
of cosmic baryon density over which the observed abun- 
dances of the light isotopes, D, 3 He, 4 He, and 7 Li agree 
with their calculated abundances. This narrow range 
in the theory's single free parameter (once the input 
physics is specified) was summarized in 1995, in units 
of critical density, as n B h 2 = .009-.020 (|. (h is the 
Hubble constant in units of 100 km/s/Mpc; this limit is 
also customarily quoted in terms of the baryon-to-photon 
number ratio, 2.5 x 10~ 10 < r] < 6 x 10 -10 .) 

The situation has changed dramatically since that 
time, with the arrival of more precise astronomical 
measurements of D, He, and Li abundances. Of par- 
ticular note are the precise measurements of the deu- 
terium abundances in high-redshift quasar absorption 
systems |^,|| . While most of the deuterium in the solar 
neighborhood has been subject to destruction in pre- 
main-sequence stars, the composition of these objects is 
believed to be nearly primordial. (Deuterium is so weakly 
bound that BBN is the only realistic site for cosmic 
production Q.) Because the amount of deuterium pro- 
duced in BBN depends strongly on baryon density, these 
measurements allow a tight constraint on that parameter 
— presently at a level of 8%, based on the measurements 
of Buries and Tytler ^,|) and the standard estimation 
of theoretical errors. In a few years, the deuterium- 
inferred baryon density will be subject to comparison 
with a similarly precise inference of the baryon density 
from observations of the cosmic background radiation || . 
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FIG. 1. The reaction network that determines yields in 
standard BBN. 



This is an important comparison, because the physical 
bases of these two inferences are completely independent. 

The new state of affairs (with corresponding ad- 
vances in He and Li observations) has been described 
as a "precision era" for big-bang nucleosynthesis [[|. 
However, as the observational uncertainties shrink, the 
uncertainties on the calculated abundances begin to 
dominate. Although the nucleosynthesis calculation itself 
is straightforward, and has been well-understood for 
over three decades now p|-|ll| , "theoretical" uncertainties 
arise from its nuclear cross section inputs. These inputs 
consist of cross sections (Fig. 1) which have been 
measured in nuclear laboratories since the 1930's, with 
accurate enough data for nucleosynthesis work dating 
mostly from the 1950's and 1960's. It is well-known that 
stellar nucleosynthesis requires cross sections below the 
Coulomb barrier, which cannot be directly probed in the 
laboratory because the cross sections are too small. In 
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contrast, BBN occurs at sufficiently high temperatures 
(T ~ f 9 K) that this is not a serious problem, and data 
are generally present at exactly the energies where they 
are needed for a BBN abundance calculation without 
any recourse to theoretical modeling. In fact, there is a 
sufficiently large body of precise data in the energy range 
needed for BBN that the results of the calculation have 
remained nearly the same since the standard code was 
first run by Wagoner in 1967 Q, despite a slow trickle 
of new cross section measurements. Nonetheless, some 
uncertainty remains in the calculations, arising from 
experimental uncertainties in determining these cross 
sections. These uncertainties range from about 5% to 
25% in the cross sections, and propagate to errors of up 
to nearly a factor of two (from lower to upper 2a limits) 
in the case of the calculated 7 Li abundance. (This is 
true of "standard BBN," which contains nothing beyond 
the standard model of particle physics. There are, of 
course, greater uncertainties if alternative - e.g. baryon- 
inhomogeneous - scenarios are considered Jl2| , ^3[ .) 

The present "industry-standard" error estimation for 
the inputs was done by Smith, Kawano, and Malaney 
(hereafter SKM, Ref. Q) in 1993, using the Monte Carlo 
error propagation method applied earlier by Krauss and 
Romanelli M). This was a landmark work because it 
examined all the nuclear inputs critically and in detail, 
and it placed quantitative error estimates on all the 
inputs. There has been a small number of new cross 
section measurements since p^Jl7[|, which have been 
unevenly incorporated in subsequent work. Subsequent 
authors have continued to use the SKM rates and errors 
for most or all reactions. A few have substituted the 
results of a single new measurement in place of the 
corresponding SKM evaluation of all experiments for the 
given reaction. 

Because the SKM uncertainties are used so widely to 
draw quantitative conclusions concerning many aspects 
of cosmology and particle physics, it is important to 
examine their assumptions and attempt improvement as 
well as to maintain some up-to-date set of uncertainties. 
The SKM work proceeds as follows: Cross section data 
for the key reactions (Fig. [I]) are gathered from an 
extensive survey of the literature. These data are fitted to 
standard, if very approximate, theoretical expectations 
concerning their energy dependence (typically a low- 
order polynomial, for the low-energy S- factor), which 
has been widely used in previous work p8|-pl|. Some 
arbitrary choices are made concerning inclusion and 
exclusion of data points in the fits, with a view to making 
these fits accurate over the energy range critical for BBN. 
Error estimation begins with formal error estimates on 
the fitted parameters, but is modified where necessary 
to include most or all data points within two-sigma 
error curves. This conservative approach justifies some 
arbitrariness in weighting the data sets. The SKM 
errors on the cross sections are summarized (with two 
exceptions) as uncertainties in the overall normalization, 
and then propagated through the basic BBN calculation 



by a Monte Carlo procedure with normalization values 
drawn from Gaussian distributions. 

We begin by listing principles for a better treatment 
of nuclear data in a BBN calculation, commensurate 
with the goals of ensuring that the confidence limits 
being quoted in cosmological work are meaningful, and 
of making a direct link to the nuclear measurements. 
Because the energy dependences of the cross sections are 
not in question in most cases, our main goal is to arrive 
at a suitable method of estimating errors in the BBN 
yields. The desirable qualities for such a method are 
are: 1) Nuclear data, and their published uncertainties, 
should be incorporated into the BBN calculation as 
directly as possible, so that errors in the calculation are 
directly linked to the nuclear data set, and the results of 
the BBN calculation reflect as few arbitrary choices as 
possible; 2) minimal assumptions concerning functional 
forms of cross sections (as functions of energy) should 
be used, where there are enough data to characterize 
these functional forms (almost all cases); 3) there should 
be accounting for correlated errors in the data, since 
they are ubiquitous; 4) the incorporation of future data 
with smaller formal errors into the inputs should reduce 
the uncertainty estimates in the calculation; and 5) 
the prescription should be no more conservative than 
necessary, given the precision with which abundances 
are now being measured. The procedure of SKM, 
outlined above, does reasonably well on point number 
2. Point number 5 is something of a matter of taste 
(although there is a mismatch in levels of conservatism 
between SKM and the errors quoted by astronomers). 
However, the SKM procedure does not answer well to 
our other requirements. The subjective assignment of 
errors causes trouble on points 1, 3, and 4. In particular, 
the conservative SKM method is not re-applicable in 
a way that narrows error estimates as new data are 
incorporated in the data set, unless old data are thrown 
out, and so it does not satisfy point number 4. Poor 
performance on this point discourages the production of 
more precise measurements, because it is unclear how 
new data will affect error estimates on reaction rates. 

To improve on the SKM work, we propose a new 
prescription for treating the nuclear data that uses the 
measured cross sections directly in the standard big- 
bang nucleosynthesis calculation. It is a Monte Carlo 
technique that uses "realizations" of the full nuclear data 
set for the key reactions identified by SKM and Krauss 
and Romanelli to derive reaction rates. We extract 
final nucleosynthesis yields and confidence limits from 
distributions of yields corresponding to the distributions 
of cross sections implied in quoted experimental errors. 
While there may be some question of uniqueness, our 
method does provide a useful and unambiguous prescrip- 
tion for linking the results of the BBN calculation to 
its nuclear inputs. This method is very much in accord 
with the stated goals of earlier work, both of SKM and 
of Krauss and Romanelli jlO| , and it reflects increasing 
computer speed. In the words of Krauss and Romanelli, 



2 



"unless error estimates are tied to a direct analysis of 
the data one cannot accurately gauge the margin for 
improvement." 

The rest of this paper is organized as follows: In 
Section |J we describe our Monte Carlo method for 
computing light-element yields and confidence intervals 
from the nuclear data. In Section III, we discuss the 



nuclear inputs in detail, emphasizing differences from the 
analysis of SKM. In Section IV, we discuss the results of a 



standard BBN calculation using our method. In Section 
we discuss the implications of our work both for the 
BBN calculation and for cosmology in general. 



II. METHOD 

For our purposes, the nuclear inputs for a BBN 
calculation come in the form of angle-integrated cross 
sections, &(E), where E is reaction energy, reported in 
the published literature. The quantities needed to evolve 
abundances in a reaction network are thermally-averaged 
cross sections, 



TTfl(kT) 3 



J a(E)Ee- E / kT dE, 



(1) 



where /i is the reduced mass and k is Boltzmann's 
constant. 

In the case of charged particles, the procedure typically 
followed in the past to obtain these thermally-averaged 
rates has been as follows: One first re-parameterizes the 
measured cross sections as S'-factors, 



S(E) ^Ea(E)e 2 ^, 



(2) 



where £ = Z\Z2a/v 1 to remove the strong Coulomb 
dependence of the cross section. Here, a is the fine 
structure constant, Zi are atomic numbers of the nuclei, 
and v is their relative velocity. The resulting function 
is then fitted to a low-order polynomial plus resonant 
terms to obtain a form which can be integrated. (This 
functional form has sometimes been used for extrapola- 
tion to low energy, although that is discouraged.) The 
integration of the nonresonant rate is then performed 
using a saddle-point approximation with lowest-order 
corrections Jl8| , |l9| ], or else performed numerically and 
fitted to the functional form one gets from the saddle- 
point integration. (The same considerations apply for 
neutron-induced cross sections, with the exception that 
they are fitted by a low-order polynomial in velocity 
before integration, and the integrations are then exact.) 
Resonant cross sections are fitted to Breit-Wigner or 
single-level i?-matrix forms, integrated numerically, and 
fitted to analytic forms. While these functional forms 
are useful for disseminating evaluated reaction rates 
in printed form, and have therefore become somewhat 
standard, they are not ideal for producing rates valid 
over a wide range of temperatures or for arbitrary S(E). 



Both SKM and Krauss and Romanelli used this sort 
of procedure, including most of the extant data in the 
fits, and being careful that the fits were performed over 
the energy range needed for BBN calculations. They 
followed it up by estimating uncertainties on the reaction 
rates, expressed as errors in overall normalization (except 
for two cases of energy-dependent errors in SKM) , which 
are appropriate to the BBN energy range. They then 
estimated the corresponding uncertainties in BBN yields 
by varying the rates according to these error estimates 
in Monte Carlo BBN yield calculations and examining 
the distribution of output abundances. The Monte Carlo 
approach was originally deemed necessary because it is 
not clear without doing the calculation whether the errors 
combine linearly. The recent work of Fiorentini et al. 
p2j indicates that in fact normalization errors on the 
reaction rates can be propagated linearly through the 
BBN calculation, and the results are very close to the 
corresponding Monte Carlo results. 

Our method differs from that of these earlier efforts 
in several ways: we fold the process of characterizing 
reaction rates into the same Monte Carlo process that 
calculates abundances, we try to make the least possible 
number of assumptions about the functional forms of 
cross-section energy dependences, and we try to keep 
error estimation based strictly on quoted errors in the 
nuclear data set, including correlated errors explicitly. 

Our Monte Carlo calculations sampled the space of 
nuclear cross sections indicated by the nuclear data, and 
the result was a distribution of output abundances, from 
which we derived confidence limits. Each calculation 
sampled 25,000 points in this space. (This was a number 
that gave us smooth 95% cl curves as a function of baryon 
density.) A single Monte Carlo sample consisted of the 
following steps: 

1. For every measured cross section of every reaction 
(every a(E) in our database), a random number 
was drawn from a Gaussian distribution whose 
mean is the reported cross section at that energy, 
and whose variance is the reported variance for 
that point. Also, for each data set (collections of 
cr(E) values from a given experiment), a random 
number was drawn from a Gaussian distribution 
whose variance was the normalization error shared 
by those points, and all cross sections in that data 
set were multiplied by this random normalization. 
This provided accounting for correlated errors. 
The database of cross sections and uncertainties is 



described in Sec. Ill below. 



2. After "synthetic" cross section data were chosen, 
smooth representations of these data were cre- 
ated for integration. Specifically, the S(E) (or 
<jv, for neutron- induced processes) curve for each 
reaction was fitted to a piecewise polynomial (in 
£>-spline representation) as a function of energy. 
The spline broke up the energy axis into several 
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(typically less than ten) segments, generally evenly- 
spaced in logE, assumed a polynomial of order 3- 
5 within each segment, and forced continuity of 
derivatives across the segment boundaries. This 
curve was fitted to the simulated data by the 
customary weighted linear least-squares technique. 
Note that this approach is "theory-free" in the 
sense that the only assumption we have made 
about the cross sections is that they are sufficiently 
smooth functions of energy to be represented by 
the chosen piecewise polynomial. The S factor 
re-parameterization is only a re-parameterization, 
and is converted back to cross sections after splines 
are fitted. Although the number of variables for 
our smooth representation was chosen by eye for 
each reaction individually, we expect that small 
variations in the functional form of the S factor 
which depend on choices made in the fitting are 
smoothed out by the subsequent thermal averaging 
and Monte Carlo sampling. 

3. The smooth representation of each cross section 
was integrated numerically at many different tem- 
peratures (ten per decade), and the resulting re- 
action rates fitted to a smooth representation as a 
function of temperature. This allowed subsequent 
calculation of rates by interpolation without expen- 
sive integrations. 

4. Finally, these reaction rates were fed into a stan- 
dard BBN code |^,[| , which used them to calculate 
yields. Note that since independent integrations 
of different numbers are performed for each Monte 
Carlo step, any (non-systematic) integration errors 
are accounted for in the Monte Carlo process. 

When the code was finished, we extracted 95% confi- 
dence limits from the distribution of yields. This process 
is similar to a treatment of fitting errors described in 
Numerical Recipes [£3| . 

We treated three reactions somewhat differently from 
the rest. First, we used the most recent experimental 
value for the neutron lifetime, 885. 4± 2 s, to derive values 
for the neutron-proton interconversion rates, choosing 
Monte Carlo values for this rate according to a Gaussian 
distribution with this mean and variance. These rates 
affect only the final 4 He abundance significantly. This is 
also the only process 4 He yields are sensitive to at likely 
values of baryon density. Lopez and Turner p4j] have 
incorporated in a coherent and consistent way a number 
of small but important physics and numerical corrections 
to the weak rates, and we take yields for this nuclide 
from their work. Second, there was not enough data 
coverage to use our technique for proton-neutron capture, 
p(n, "f)d, so we used a theoretical model of this process 
as described in Sec. 
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below. Finally, the apparent 
presence of systematic discrepancies in measurements 
of the 3 He(a,7) 7 Be cross section required the special 



treatment described in Sec. [II D 8, after applying our 
standard technique to a subset of the data. 



III. NUCLEAR INPUTS 
A. The Database 

The nuclear data used in our work were obtained from 
a comprehensive survey of the experimental literature 
from approximately 1945 onward. Many of the numerical 
values were obtained from the on-line CSISRS database 
p5[ . However, even in these cases, we incorporated data 
sets only after reading the original sources carefully. For 
almost all reactions, the data sets included were the 
same ones found in SKM. There are three reactions, 
d(p, 7) 3 He, 3 He(n,p) 3 H, and t(a, j) 7 Li, for which more 
recent cross section data than those used by SKM exist. 
Some subsequent BBN calculations (e.g., Ref. have 
incorporated the latest measurement of t(a,^f) 7 Li |1J] by 
replacing the SKM fit with a reaction rate and uncer- 
tainty based on that measurement alone. These numbers 
have not found their way into all subsequent work (e.g., 
Ref. ^2|), and we know of no calculations that have 
incorporated cither of the recent TUNL measurements 
of the d(p, 7) 3 He |jl5],[l6) cross section (save one reported 
in this last experimental reference). One of our goals is 
to provide a new standard calculation which incorporates 
these measurements so that they are not ignored in future 
theoretical work. 

Another of our major goals is to incorporate explicitly 
information concerning known systematic errors in the 
cross section measurements, which indicate correlations 
among data from a given experiment. The incorporation 
of this information into our Monte Carlo calculation 
is described in Sec. [n] above. Here, we describe the 
origin of the numbers used. Data sets which include 
measurements of the cross section for the same process at 
several energies contain both shared normalization errors 
(which affect all points in that data set) and unshared 
point-to-point errors (which vary from one point to 
the next, and arise primarily from counting statistics). 
Because it is generally much more difficult to determine 
the absolute normalization of a cross section than it is to 
measure its energy dependence, most of the uncertainty 
in a given cross section is usually in its normalization. 
It is therefore important to treat normalization errors in 
the data properly when combining data sets, especially 
if the data sets are of different sizes. It is not correct 
to apply the quoted "total error" for each point in a 
fit when the correlated errors have been added in (the 
usual case for published data). Our calculation assumes 
Gaussian-distributed errors as the simplest assumption in 
every case, although one might expect that normalization 
errors, in particular, will not be Gaussian-distributed. 
Note that in other contexts, one typically allows a 
floating overall normalization for each data set in fitting 
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nuclear data. Since the energy dependences are well- 
determined for almost all of the eleven key cross sections, 
we are only interested in experiments which measure the 
cross section normalization. 

Separate treatment of shared and unshared errors 
involves some difficulties. These include the fact that 
systematic effects may not be well-quantified, and the 
fact that experimental data have often been presented 
in ways incompatible with this goal. Several of the 
experimental sources (especially the more recent ones) 
explicitly state shared normalization errors separately 
from unshared errors. In some cases, quoted normal- 
ization errors had to be subtracted from quoted total 
errors, which are quadrature sums with unshared errors. 
Other experimental sources did not explicitly add up 
the shared errors, but they usually listed the percent 
contribution from each major item in the error budget 
separately. These could usually be added up to provide 
normalization errors, with a small amount of guesswork 
as to which contributions to place in each category. 
For example, detector efficiencies and target chemistry 
are often readily identifiable as sources of shared error, 
in those cases where they are. Because errors are 
almost always added in quadrature, our procedure was 
not strongly dependent on identifying which category 
(shared or unshared) each individual contribution should 
be placed in, so long as we were careful that estimates 
of unshared errors did not become too small. We are 
confident that the sizes of errors adopted for our database 
are all roughly correct. Our method required that we 
exclude a small number of data sets which did not 
provide enough information for such a breakdown of error 
sources (but we were biased toward keeping data sets, 
especially where there were few measurements of a cross 
section). We also excluded data sets which only represent 
measurements of relative cross sections. Fortunately, we 
did not need to exclude a large number of data sets for 
any reason. 

We have not incorporated any of the substantial body 
of theoretical knowledge (ranging from unitarity con- 
straints to "microscopic" reaction models) which exists 
for some of these processes into the database, except for 
the process p(n,j)d. The reason for this is that there 
is little need for theoretical evaluations where we already 
have large amounts of data, especially since errors are of- 
ten hard to assign to models. A typical use of theoretical 
models is to calculate an energy dependence for a cross 
section, test it by comparison with relative cross section 
measurements, and then normalize it by absolute cross 
section measurements. This is more constraining than 
our approach, since it allows data taken at all energy 
ranges (and sometimes in other reaction and scattering 
channels) to affect the evaluated cross section at a given 
energy, but it would require a significantly larger and less 
straightforward effort. In this regard, our "theory-free" 
approach may still provide conservative error estimates 
for some reactions. 

Because of all these concerns, particularly those con- 



cerning unknown systematics and arbitrariness of error 
assignment in experiments, there can probably be no 
unique evaluation of rates and their errors. We have 
attempted to make estimates with the desirable qualities 
listed above, and to improve on previous efforts. 

B. Reaction Sensitivities 

Because our approach to the BBN calculation is very 
closely tied to the data, we are able to make very precise 
statements about exactly which inputs the calculated 
yields are sensitive to (at least for the "standard BBN" 
calculation we have done). In particular, we have 
quantified this sensitivity with what we call "sensitivity 
functions" for each reaction and nuclide. These may be 
regarded as functional derivatives of BBN yields with 
respect to reaction S factors. For a given reaction and 
nuclide, we calculate one of these functions numerically 
by adding a small amount to the reaction S factor over 
a narrow bin in energy, and computing the resulting 
primordial abundance of that nuclide. The sensitivity 
function is the fractional difference between this yield 
and the unperturbed yield, as a function of the location 
of the energy bin. In the limit of very narrow energy bins, 
this is a functional derivative. The functions indicate 
quantitatively the exact energies at which each process 
is important in standard BBN, and therefore the exact 
energies at which precise cross section measurements are 
needed to produce precise BBN calculations. We denote 
the sensitivity functions by g2(E) for D/H and gi{E) for 
7 Li/H, and they are shown for all eleven cross sections 
below. 

The general behavior of these functions is clear: Above 
some temperature in the vicinity of 10 9 K, the actual 
reaction rates do not matter because the reactions are in 
thermal equilibrium, with forward and reverse reactions 
proceeding at the same rate. When the density and 
temperature drop to some point at which a particular 
reaction falls out of equilibrium, its rate may become 
a determinant of the nuclide abundances. Finally, 
there comes a point when temperatures and densities 
are too low for a reaction to change abundances at 
all. The approximate effective energies of reactions for 
nuclear burning at a given temperature correspond to 
the customary "Gamow peak" of nuclear astrophysics 
]l8| , and what we see in the sensitivity function is often 
the Gamow peak, convolved with the distribution of 
temperatures at which a given reaction is active, but out 
of equilibrium. The rate sensitivities are also functions 
of th e baryon density, as indicated in the plots of Sec. 

hid. 
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C. Statistical tests 

Our new method for calculating BBN abundances 
should be examined to test the robustness of its results 
and the consistency of the input nuclear data. We 
therefore applied it to fake data drawn from the assumed 
source distribution of the actual data. We gener- 
ated (for each reaction) fake data that were Gaussian- 
distributed about our highest-probability curve. Cor- 
related data were multiplied by appropriate Gaussian- 
distributed common normalizations. The fake data were 
placed at the same energies as the actual data, and 
the distributions were based on the quoted errors of 
the corresponding actual data. For each collection of 
fake data, we applied our Monte Carlo curve-fitting 
method, and computed the means, £t(6j), and standard 
deviations, a(bi), of the _B-splinc coefficients describing 
the fitted curves, (bi denotes the coefficient of the ith 
basis function. Note that the local support of B-spline 
basis functions [^6| makes each of these coefficients a 
sort of weighted local average of the fitted function, 
so that fJ>(bi) and o~(bi) reflect the means and standard 
deviations of the fitted curves directly.) After performing 
this procedure for many (> 1000) collections of fake 
data, we computed the means, o~(bi) and standard 

deviations E(cr(6j)) over the fake- data distribu- 

tion of the -B-spline coefficient means and standard 
deviations o~{bi) derived from each choice of fake data. 

Comparison of these numbers with the means and 
variances of -B-spline coefficients generated as interme- 
diate numbers by our BBN code indicates that the 
means reproduce the assumed curves consistently with 
our error estimates in all but one case. The mean 
variances a(bi) of the B-spline coefficients generally agree 
with the standard deviations of .B-spline coefficients from 
our BBN procedure (which are o~i for the actual data); 
where they differ, our BBN procedure gives slightly larger 
standard deviations. The standard deviations of the 
means, were virtually always identical to the 

mean standard deviations, a(bi). In other words, the 
distributions of S-factor curves from our procedure do 
not change drastically when the experimental data are 
drawn from other points in the same distributions from 
which we assume the actual data to be drawn. Serious 
errors in reproducing the assumed curves occurred only 
for 3 He(d,p) 4 He. 

We repeated the same numerical experiment, this 
time simulating a factor-of-two underestimation of all 
normalization errors by adding extra scatter to the fake 
data (doubling all of the normalization errors). In 
applying our Monte Carlo method to these fake data, 
we assumed the quoted normalization errors, not the 
inflated ones. The most noticeable results were an 
increase in the standard deviations of means S(/i(6j)) by 
a factor of about y/2, and slighly poorer reconstruction 
of several coefficients for t{a,"f) 7 \A. Very little else 
changed. The standard deviations over fake data of 



the -B-spline variances, S(cr(6,)), also increased in some 
cases, but were seldom much more than 10% of cr(bi). In 
other words, the variances of the S-factor curves from 
our BBN procedure are sensitive mainly to the error 
estimates rather than the quoted cross sections, while 
the mean curves are less sensitive to the error estimates. 
Underestimated errors result in S'-factor curves that 
vary less than they should in our BBN calculation, but 
probably never by factors of more than 1.5. 

Although the results above suggest that our method 
is not overly sensitive to unexpected scatter in the data, 
we still wish to address the question of consistency of 
the nuclear data. We calculate a chi-squared statistic 
with the data and the "best-fit" model curve describing 
the un-altered data. This is not a goodness-of-fit test 
for our modelling, because the individual curve plays 
no important role in our calculation. It is, rather, a 
benchmark by which to examine the consistency of data 
measured at different energies. Initially, we calculated 
this chi-squared statistic using the full non-diagonal error 
matrix, as described in Ref. |2^ ]. However, as this 
reference shows, chi-squared calculations using the non- 
diagonal error matrix are not appropriate when the 
correlations are in the form of shared normalizations. 
Fits minimizing this statistic are almost always lower 
than the data when there are large normalization errors. 

The statistic that we did apply is the ordinary \ 2 
statistic, defined as 

2 _ V"^ (S(Ej) - ffmodcl^)) 2 ,„s 

X — a 2 > \ 6 > 

i 1 

where S(Ei) are the data for a given reaction, measured 
at energies Ei, S mo dci(E) is the best-fit model curve, o~i is 
the quadrature sum of normalization and point-to-point 
errors for each datum, and the sum is over all data for 
a given reaction. Because of the correlated errors, we 
do not expect this statistic to be distributed as a formal 
chi-squared distribution. We therefore constructed the 
expected distribution by varying fake data about the 
best-fit curve according to the quoted uncertainties as 
above. This allowed us to assess consistency of the data 
by determining the likelihood of the actual value of % 2 > 
given our assumptions about the distribution from which 
the actual data were drawn. The results of computing 
X 2 for each reaction S'-factor and comparing it to the 
fake data distribution are shown in Table ffl. From these 
results, we conclude that the only clear case of trouble 
is 3 He(d,p) 4 He. Note that this cross section has no 
noticeble effect on the calculated deuterium abundance, 
and very little effect on the calculated 7 Li abundance. It 
is more important for 3 He, which is not as yet amenable 
to high-precision treatment as a product of BBN. 

The reader may object that we do not have a procedure 
to inflate the errors where the \ 2 test indicates that 
they may have been underestimated, or to discard data 
sets which contribute disproportionately to x 2 ■ We have 
resisted adopting such procedures for several reasons. 
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TABLE I. For each reaction, the number of data points, characteristics of the piecewise polynomial descriptions, values of 
\ 2 , and probabilities that \ 2 exceeds its measured values (based on 20,000 Monte Carlo samples). 







Spline 






Probability 


Reaction 


Points 


segments 


order 


2 

X 


of exceeding \ 


p(n, 7 )d 




1 


10 




(not fitted to lab data) 


d(p, 7 ) 3 He a 


16 


5 


3 


3.40 


93% 


d(d,p) 3 H 


134 


6 


3 


108.0 


68.7% 


d(d,n) 3 He 


130 


6 


3 


106.7 


64% 


3 He(a,7) 7 Be b 
3 He(d,p) 4 He 


118 


4 


3 


138.3 


37% 


111 


9 


4 


315.0 


0.5% 


3 He(n,p) 3 H 
7 Li(p,a) 4 He 


240 


5 


5 


106.0 


91% 


105 


6 


3 


193.1 


11% 


7 Li(p,n) 7 Be 


137 


29 


3 


223.0 


16% 


3 H(a, 7 ) 7 Li 
3 H(d,n) 4 He 


55 


6 


3 


129.5 


12% 


213 


15 


3 


292.8 


13% 



a This reaction is a special case, since none of the experiments in our database overlap in energy. 
b Only capture photon data are described here. 



The first is that it would be difficult to arrive at a 
unique and meaningful prescription to determine by how 
much to inflate errors, or which data to discard. (The 
customary multiplication by the square root of the x 2 
per degree of freedom is not suitable because \ 2 IS n °t 
distributed as a formal x 2 distribution, and because 
the factor should in principle be a function of energy.) 
The second reason is that we would like to leave any 
systematic problems with the existing data, and their 
properties with regard to generating formal errors for 
BBN, as distinct issues. A database which performs 
well on both counts is desirable for producing reliable 
BBN calculations. On the other hand, discarding data 
is something we will not do without good reasons from 
the nuclear physics literature, even though disproportion- 
ately large fractions of the \ 2 statistics usually do come 
from individual data sets in our database. Finally, there 
is only one "two-sigma" unlikely x 2 m Table [| although 
there are more "one-sigma" unlikely values of % 2 than 
expected. 



D. Individual Reactions 

The following section is dedicated to detailed dis- 
cussions of the data sets for each individual reaction. 
For convenience, we have followed the same ordering 
as SKM in discussing individual reactions. Because the 
available data have changed very little since the work of 
SKM, we emphasize differences from their analysis and 
new insights provided by our examination of the inputs. 
We also list the references from which we compiled 
our database. Unless explicitly discussed, any omission 
was because either a normalization error could not be 
extracted from the original source, or the experiment 
measured only relative cross sections (or, trivially, the 
experiment contained no data in the relevant energy 
range) . 



For each reaction, we show in a graph (Figs. |^-|l^) 
the input data, the "best-fit" curves and 95% confidence 
limits inferred from our Monte Carlo procedure (solid 
curves), and (where available) the corresponding curves 
from the SKM analysis (dash-dotted curves), and curves 
from the ENDF-B/VI [£8| evaluation (dashed curves). 
The tick marks below the data in each graph show the 
region over which the integral must be performed to get 
the yields correct to one-tenth of the total uncertainty in 
all abundances (inner tick marks), and to get the yields 
correct to one part in 10 5 (outer tick marks). 

In the lower panels of each graph, we show the 
sensitivity functions at baryon densities of il^/i 2 = 0.019 
(solid curves) and Vtsh 2 — 0.009 (dashed curves) for D 
and 7 Li. They represent very well the contribution of 
the S factor uncertainties to the final yield uncertainties 
as functions of reaction energy in the sense that, to a 
good approximation, the convolutions of these curves 
with the 95% cross section limits gives the 95% yield 
limits due to that reaction indicated by the Monte Carlo 
study described below. The g functions shown in the 
plots have been multiplied by energy in MeV so that 
relative areas under the curves can be accurately judged 
on the logarithmic scale on which we have plotted the 
nuclear data. Some references contained two or more 
distinct data sets with different normalization errors; 
in these cases, an extra number after publication year 
distinguishes symbols from the same publication. All 
energies are in the center-of-mass frame. 



1. p(n,f)d 

As noted above, this reaction required special treat- 
ment. It could not be subjected to our piecewise splining 
or Monte Carlo sampling of data points because of an 
extreme scarcity of data in the Ecu < 300 keV range. 
This is despite recent efforts to measure this crucial cross 
section [B0By| in the relevant energy range, which have 
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FIG. 2. Theoretical curve from Hale et al. |2J, SKM fit to 
that curve with 2a errors, and current experimental data in 
the BBN energy range J3cj , [jl| . An explanation of our nuclear 
data figures may be found near the beginning of Sec. [II D 



so far resulted in only four data points. We have chosen 
as the cross section for this process its evaluation in the 
ENDF-B/VI database |29|. This cross section derives 
from a theoretical model computed sometime around 
1970, constrained by np scattering phase shifts, along 
with a few photodisintegration data (known to have 
systematic problems), and it has seen a minor update to 
match the current value of the well-measured thermal- 
neutron capture cross section No documentation 
for this model survives. The crucial energy range for 
BBN corresponds to a changeover in reaction mechanisms 
(from Ml capture to El capture) for this process, so one 
might expect the validity of the evaluation to be most in 
question there. It has held up remarkably well in light 
of the recent measurements, but its authors view this 
as a fortuitous coincidence |33|. SKM estimated a 5% 
uncertainty on this evaluation, based on uncertainties 
quoted in earlier tabulations of the cross section for 
practical use. It is very difficult to trace the origin 
of this number, but we have adopted a 5% Gaussian- 
distributed normalization error as the uncertainty in 
this cross section — both for consistency with SKM, 
and because it is not too far from an estimate by the 
evaluation's authors of "at least 10%" Q. (Note that 
the 7% total uncertainty quoted by SKM includes errors 
in further fitting and integration. We do no further 
fitting, and we estimate an error of less than 1% in our 
numerical rate integrations.) 



2. d(p,j) 3 Re 

The data set for this reaction is not well suited to 
our method. Of the three experiments performed before 
1990, the experiment of Bailey et al. |54| is only a relative 
measurement, and not suited to our method; the low- 
energy experiment of Griffiths et al. p5| deconvolved 
thick-target data with a simple model of combined S- 
and P-wave capture at a sharp nuclear surface. The 
resulting energy dependence is in conflict with both the 
more recent measurements of Schmid et al. JT^ j and the 
three-body microscopic calculation of Viviani et al. |56| . 
The recent TUNL experiments |H||l6| also suggest that 
these earlier thick ice target experiments used the wrong 
stopping powers, resulting in cross sections about 15% 
too high. We excluded the Griffiths et al. and Bailey et 
al. data sets for these reasons, but only after checking 
that their omission did not alter our results drastically. 
As indicated below, this reaction contributes a large 
portion of our uncertainty estimates. The only other 
experiment in our database for this reaction is from Rcf. 



The sparseness of the data may be a cause for concern 
regarding our treatment of this cross section. However, 
it is encouraging that our splines indicate much the same 
energy dependence as the microscopic calculation |5(||l6| , 
so the sparseness is probably not such a big problem. Our 



8 




errors reflect the ~ 8% normalization errors in the best 
measurements of the cross section. 



3. d(d,n) 3 Re 

The deuteron-deuteron reactions have been measured 
very extensively for fusion applications, and they are 
especially well-constrained below about 60 keV by the 
precise measurements of Brown and Jarmie jB8|. The 
SKM errors were based on an apparent discrepancy 
between these data and those of Krauss et al. . Brown 
and Jarmie make a similar assessment of the situation, 
and recommend renormalizing the Krauss et al. data 
to their more precise data (which we do not do). It is 
important to note that most of the error in both cases 
is contained in normalizations, so the case is not a 10% 
discrepancy between ten data points with 8% errors and 
eleven points with 1.5% errors. Note that the sensitivity 
functions peak at 100 keV or above, beyond the Brown 
and Jarmie data. With the exception of Krauss et al. 
(with an 8% claimed normalization error), all of the data 
in this region date from before 1960, and their large errors 
and scatter are responsible for the large contribution of 
this reaction to the uncertainties in BBN yields. Other 
data for this reaction are from Refs. |4C 



4- d(d,p)t 

This cross section is generally measured concurrently 
with that of d(d, n) 3 He, so most of the discussion for that 
reaction carries over to this one. The sensitivity functions 
here are more complicated, but there is still a large 
contribution at energies greater than 100 keV. There are 
also fewer data above 100 keV in this case, reflecting 
the use of neutron-specific detection methods in some 
d(d, n) 3 He experiments. Again, few of the data which 
are present in this range come from modern experiments. 
The data for this reaction were taken from Refs. fB8r-H. 



5. 3 He(n,p)i 

The data set for this reaction consists of cross sections 
for both the forward (Refs. pE| - |48| ) and reverse (Ref. 

processes. We used reverse data that had 
been converted to forward cross sections through exact 
detailed balance relations, along with direct data. The 
SKM fit included numbers from Alfimenkov |50|, which 
we did not obtain, but the important difference between 
our analysis and that of SKM is that they excluded the 
very precise measurement of Borzakov et al. p5| , which 
owes its small quoted uncertainty to a normalization from 
6 Li(n, a) 3 He data intended for metrological use |pl|,p2|. 
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FIG. 5. S factor data and fits for d{d,p)t. 
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FIG. 6. S factor data and fits for 3 He(n,p)t. 



SKM excluded this data set because it decreases with 
energy more quickly than the other data, ending lower 
than any of them at 100 keV. This, combined with the 
small quoted errors, would have forced their second- 
order polynomial fit below the data at higher energies. 
However, this is not a problem for our piecewise spline 
approach, which decouples errors arising from different 
experiments at different energies. The extensive and 
much-needed data which Brune et al. [ft7| published 
recently for this cross section also agree more closely with 
the Borzakov et al. than with SKM. 



6. i(d,n) 4 He 

The definitive measurement of the cross section for 
this process, up to 70 keV, is by the Los Alamos group 
f^pi! , and it has a quoted normalization error of 1.4%. 
Other measurements in our database come from Refs. 

As SKM point out, there are no modern 



mi 

experiments on the high-energy side of the resonance. 
However, this reaction does not contribute noticeably to 
the uncertainty in D or 7 Li yields. 



7. 3 He(d,p) 4 He 

There is considerable scatter in the data for this 
process, as reflected in the very low probability of the \ 2 
statistic. In the case of the Bonner et al. measurement 
p9[ , this may be attributable to energy straggling in the 
gas target's Al window. In any case, there is little else 
to say here except to point out that our piecewise spline 
seems to represent the data about as well as the SKM 
i?-matrix fit. This reaction contributes a large portion 
of the uncertainty in the 3 He yield. Other data used for 
this reaction are from Refs. |3£j and J60|-|62[. The low- 
energy data of Schroder et al. |pq| , intended to probe 
electron screening, were omitted here because they did 
not measure the absolute cross section. 



8. 3 He(a,7) 7 Be 

The apparent discrepancy at 2a between experiments 
that observe capture gamma rays from this process (Refs. 

[j9|) and those that observe photons from the 7 Be 
decay (Refs. p^-fnl) is well-known from discussions of 
the solar neutrino problem j72|. Our direct approach 
to the data makes no use of the (well-understood) 
theory of this process, and the activation technique has 
only been used at energies too high to affect big-bang 
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FIG. 7. S factor data and fits for t(d, n) 4 He. 
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FIG. 8. S factor data and fits for 3 He(d,p) 4 He. 

nucleosynthesis (although the Volk et al. jnj results are 
quoted as extrapolated 5(0) values only). However, 
the energy dependence of this reaction is sufficiently 
well-determined from the capture photon experiments to 
justify the use of an "activation method" curve obtained 
by renormalizing the fit of the capture-photon-derived 
cross sections to match the activation points. 

After the main run, we re-ran our Monte Carlo sam- 
pling, renormalizing the 3 He(a,7) 7 Be cross section from 
the capture-photon measurements by the mean of the ac- 
tivation measurements, and drawing the renormalization 
from a Gaussian distribution based on the variance of 
the activation measurements. For this purpose, we use 
the reference of all cross sections to zero energy found 
in Adelberger et al. j72|. The effect of this discrepancy 
(amounting to a systematic shift of 13%) is easy to 
understand; all 7 Li produced at high comes from 
this reaction, so changes in its rate result directly in 
changes of the final 7 Li yield. The result is a shift of 
11% in our confidence limits for Li, a significant fraction 
of the widths of these limits, at high fig. (See below.) 
If the activation measurements are correct, this would 
exacerbate the problem of 7 Li depletion in halo stars — 
a point which has not previously arisen in discussions 
of this problem because it represents a much smaller 
fraction of the widths of the SKM 7 Li limits, and because 
SKM dropped the activation data from their evaluation 
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FIG. 9. S-factor data and fits for 3 He(a,7) 7 Be. Note that 
the data of Krawinkel et al. |64| have been renormalized as 
suggested by Hilgemeier et al. pS^ ■ 
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FIG. 10. S factor data and fits for i(a,7) 7 Li. 
altogether. 



9. t(a, -y) 7 Li 

A precise new measurement of the cross section for 
this process has been completed by Brune et al. |l4|] in 
the time since the publication of SKM. Some subsequent 
calculations have incorporated the reaction rate derived 
from that measurement by its authors {e.g., Copi et al. 
|TJ), using as the uncertainty the 6% uncertainty in the 
experiment's cross section normalization. Krauss and 
Romanelli pointed out that it is not always the 
best policy to base calculations on only the most recent 
measurement, and we agree. The previous data, Refs. 
]73| frq , do show a great deal of scatter, and our best- 
fit curve tends to follow the more-precise Brune et al. 
data, as it should. Relative to the SKM yields, this has 
little effect at high baryon density, but it decreases the 
7 Li yield slightly at very low baryon density, where this 
process makes most of the 7 Li. We have omitted the 
Coulomb-breakup measurement of Utsunomiya et al. |77]] 
because the Coulomb-breakup process is not completely 
understood for this reaction (as discussed in SKM) and 
because the cross section energy dependence derived from 
this method disagrees with the Brune et al. data. We 
note that the problem with normalization of the Schroder 
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et al. [Q and Griffiths et al. [79] data sets mentioned by 
SKM seems to have gone away in the intervening time 



10. 7 Be(n,p) 7 Li 

We fitted the cross section for the reverse of this pro- 
cess, since there are no direct data in the energy region of 
interest for BBN. This may not be the ideal choice, since 
the curve has a very large second derivative just above 
threshold. We fit data for 7 Li(n,p) 7 Be from Sekharan et 
al. tfr9| , Taschek and Hemmendinger | pQ| , and Gibbons 
and Macklin ]49|] , which extend from threshold to well 
past the lowest-energy resonance. We also cut off these 
data sets at about 700 keV above threshold (where they 
no longer affect BBN yields), so that only data in the 
critical range affect the x 2 calculations. Although this 
reaction contributes very little to our error budget (see 
Sec. IV below), it is important to recognize the scarcity of 
data and the lack of detailed error analysis in the original 
sources. 



11. 7 Li(p,a) 4 He 

The low-energy cross section for this process is repre- 
sented by a small number of measurements, Refs. [pl[-p5| , 
and is determined over most of the energy interval of 
interest to us by the data of Rolfs et al. Q. We 
have kept the data of Engstler et al. p4j below 50 keV 
in our fit, even though the slight rise in this data set 
relative to that of Harmon |85| has been attributed to 
electron screening effects. (The data of Harmon were 
normalized to the 6 Li(p, a) 3 He cross section in such a way 
that they have been "corrected" for screening to some 
extent [gij.) The fact that our S- factor curves follow 
the small error bars of the low-energy Engstler et al. 
data is not too distressing, since the sensitivity function 
for lithium production for this reaction is very small 
below 40 keV, and a true correction for screening would 
require a thorough theoretical treatment, e.g., i?-matrix 
techniques with a direct reaction mechanism. It is not 
clear that simply extrapolating fits to simple functions 
from higher energies is valid, or that the data of Harmon 
et al. reflect a true correction for screening. 
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E. Recommended Rates and Errors 



A disadvantage of our method, relative to earlier 
efforts, is that it requires the whole nuclear database 
and a larger amount of computer time, as well as a 
significant modification of existing BBN code. Reaction 
rate uncertainties are also harder to quote than in 
the SKM prescription, since uncertainties at different 
temperatures are neither completely correlated nor com- 
pletely uncorrelated. In any case, the rates produced 
by our code are not suitable for any use significantly 
different from the standard BBN calculation, and our 
method is specific to the BBN context, where only a few 
cross sections — well-represented at the right energies 
by available data — are needed. We note, in particular, 
that the "flaring" of our piecewise polynomial fits at high 
and low energy does not contain any physical information 
at all, but only the fact that polynomial interpolations 
always blow up beyond the limits of the data they were 
fitted to. For general use (especially outside the BBN 
energy region), we suggest using rates from a more 
general compilation, such as the NACRE compilation 
Pl , intended to succeed the Caughlin and Fowler pl| ] 
charged-particle reaction rates. 



IV. RESULTS 

In discussing the results of applying our prescription, 
we concentrate on the predictions of 7 Li and D yields. 
On the one hand, the observational status of 3 He is 
not such as to motivate precise comparisons with the 
calculation. On the other, the errors in 4 He, especially 
at higher values of f2s , are dominated by the uncertainty 
in the weak coupling constant and by uncertainties in 
calculating the matrix elements for the weak processes. 
These errors and theoretical uncertainties have been 
analyzed exhaustively by Lopez and Turner [ p4j , and we 
have not included all the apparatus of their treatment in 
our BBN code; we take their results for this nuclide to 
be definitive. 

Our most important results, apparent from Figs, [IT^, 
|l4| , and [l5[ are reductions by factors of up to three in 
the width of the 95% confidence intervals for both the 
7 Li and D yields relative to SKM. The median values of 
our yields are almost identical to those obtained from 
the SKM rates, the difference being well within our 
estimated errors. These small changes can generally not 
be attributed to any one reaction, but to some nonlinear 
addition of changes from several reactions, as indicated 
in Fig. [TJ. 

We did not expect a huge change either in the size 
of the error estimates or in the calculated most-likely 
BBN yields. We did expect some reduction in the 
error estimate because we did not go out of our way 
to be conservative. The reduction in the size of our 
error estimates is largely a side-effect of our method of 
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FIG. 13. Summary of the 95% confidence intervals for 
the BBN predictions for D, 3 He, 4 He and 7 Li. The 4 He 
uncertainty comes from Ref. E41. Boxes indicate 95% cl 
abundances from observation [p7|-[g9|] . The vertical band 
indicates our 95% cl baryon density inferred from the Buries 
and Tytler deuterium observations. 
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FIG. 14. Fractional changes (solid curves) and revised 95% 
confidence limits (shaded regions) in D /H from our full Monte 
Carlo calculation and from Monte Carlo studies of individual 
reactions, relative to the SKM yields and errors (dashed lines). 
For comparison, the vertical bar at the side of the figure 
indicates the la uncertainty in the Buries and Tytler D/H 
measurements H. All curves are on the same scale. 
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FIG. 15. Changes (solid curves) and 95% confidence in- 
tervals (shaded regions) in log( 7 Li/H) from our full Monte 
Carlo calculation and from Monte Carlo studies of individual 
reactions, relative to the SKM yields and errors (dashed 
lines). The vertical bar on the side indicates la error in the 
Spite-plateau abundance [Q. The underlined reaction label 
indicates the individual-variation result using the alternate 
rate for 3 He(a, 7) 7 Be, as described in the text. All curves are 
shown on the same scale. 
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TABLE II. Size of normalization error necessary in a 
new data set (of twenty points evenly-spaced from 5 to 
500 keV) for its addition to our database to reduce the 
single-reaction error estimate on the indicated BBN yields 
by half at Q B h 2 = 0.019. 



Reaction 


Desired error in cross section 


BBN product 


d(p,7)"He 


3.2% 


D 


'd(d,n) 3 He 


2.0% 


D 


d(d,p) 3 R 


1.6% 


D 


3 He(a,7) 7 Be 


2.3% 


7 Li 



handling the nuclear data, chosen as a way of relating 
uncertainties in the calculated yields to uncertainties in 
the nuclear data. 

We also studied the error contributions of individual 
reactions by setting all rates but one to their SKM 
values, and applying our Monte Carlo technique to that 
reaction alone. As shown in Figs. |l4| and [l5|, this 
indicates the relative importance of each reaction as a 
function of baryon density (much as in the figures of 
Krauss and Romanelli |H]]). This shows exactly which 
reaction rates need improvement to reduce the errors 
on BBN yields. In turn, the knowledge of which cross 
sections need improvement can be combined with the 
sensitivity functions of Sec. Ill and Figs. || through 
|l2| to indicate the specific energies at which they need 
improvement. The "most wanted" for the deuterium 
abundance are, from most to least important at fish 2 = 
0.019: d(d, n) 3 He above 100 keV; d(p, 7) 3 He everywhere; 
d(d,p) 3 R above 100 keV, and p(n,j)d at 30-200 keV. 
For 7 Li, the leading contributions to the uncertainty 



at Q B h 2 = 0.019 are from p(n : j)d at 20-150 keV, 
3 He(a, 7) 7 Be at 150-375 keV (and overall normalization), 
d(p, 7) 3 He everywhere, and d(d, n) 3 He above 100 keV. At 
fig/i 2 = 0.009, the uncertainty in 7 Li comes mainly from 
3 H(a, 7 ) 7 Li and 7 Li(p,a) 4 He. 

Taking this list one step further, we have generated 
fake data for some of these reactions (following the best- 
fit curve) and placed it in our database. In each case, 
fake data were placed on twenty evenly-spaced intervals 
between 5 and 500 keV center-of-mass. We then re- 
did the single-reaction Monte Carlo for that reaction, 
and reduced the size of the normalization error on the 
fake data set until the uncertainty estimate due to 
that reaction was reduced by half. The sizes of the 
normalization uncertainties required by this criterion are 
given in Table [H| We assumed that unshared errors can 
be made arbitrarily small, and left them out. Similar fake 
data sets modeling proposed experiments could be used 
to determine what effect they would have on the BBN 
error estimates in our formalism. 

While the 3 He chemical evolution and abundance mea- 
surements are too uncertain to motivate high-precision 
comparisons with the calculation §J9C]-|9§], we have also 
examined the results of our calculations for this nuclide 
and for its reaction-rate dependences. The results are 
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FIG. 16. Fractional changes (solid curves) and 95% confi- 
dence intervals (shaded regions) in 3 He/H from our full Monte 
Carlo calculation and from Monte Carlo studies of individual 
reactions, relative to the SKM yields and errors (dashed lines). 
All curves are shown on the same scale. 



indicated in Fig. [h]. The slope of the dependence of this 
abundance on baryon density has changed relative to the 
SKM rates. At low fls, this reflects the reduced rate 
of 3 He(n,p) 3 H at about 100 keV in our calculations (a 
result of our fit emphasizing more precise measurements, 
but reinforced by very recent measurements). At high 
£Ib, this reflects the reduced rate of d(p, 7) 3 He indicated 
by recent improved measurements. These effects cancel 
in the middle of the Copi et al. JD] concordance interval, 
so that the disagreement with SKM is not serious in the 
likely range for Q,b- Since much of the post-big-bang 
evolution of D/H and 3 He/H is expected to consist of 
the burning of deuterium into 3 He, the sum of these two 
number densities is often considered in comparing them 
to the BBN predictions. Therefore, we also show the 
limits on this sum from our Monte Carlo in Fig. 



V. CONCLUSIONS 

In summary, we have applied a new Monte Carlo 
approach to the use of nuclear data in big-bang nucle- 
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FIG. 17. Fractional changes (solid curves) and 95% con- 
fidence intervals (shaded regions) in ( 3 He+D)/H from our 
full Monte Carlo calculation and from Monte Carlo studies 
of individual reactions, relative to the SKM yields (dashed 
lines). All curves are shown on the same scale. 



osynthesis calculations. This approach has the virtues 
of coupling the results of BBN calculations and their 
error estimates closely to the available nuclear data, of 
explicitly handling correlated errors in the data set, of 
allowing easy use of new data, and of taking some of 
the data-evaluation process out of human hands. We 
have also abandoned the explicit conservatism of the 
previous "industry-standard" error estimation, but the 
chief virtue of our method is the close coupling of the 
BBN calculation to the nuclear data, particularly with 
regard to incorporating new data. 

Application of our method has resulted in a reduction 
in the estimated uncertainty in the BBN calculation. The 
old estimates of the uncertainties were larger than the 
quoted uncertainties in recent astronomical observations 
of BBN nuclides, so that uncertainties in the nuclear in- 
puts dominated inferences from individual observations. 
Given our prescription, the uncertainties in the calcula- 
tion are once again smaller than the quoted errors on 
any of the current observations, and this strengthens the 
constraints which can be placed on BBN. The constraints 
that can be derived by applying our calculation to current 
observations have been discussed in a previous paper fl93| . 

The most important result of our earlier paper is 
that useful inferences can now be made where nuclear 
uncertainties formerly precluded any strong conclusions. 
An example is the question of whether observed lithium 
abundances are consistent with low deuterium obser- 
vations, in the absence of lithium depletion on the 
Spite plateau. The conservative uncertainty estimates 
would also not allow any determination of the baryon 
density based on BBN to better than about 10%; our 
prescription reduces the uncertainty in this quantity so 
that even given observational errors of a few percent in 
deuterium, the uncertainty on the baryon density would 
be dominated by astronomical observations. 

We have subsequently found that the calculated BBN 
yields in Ref. J93| suffer from a programming error which 
we have now corrected. The error resulted in overes- 
timates of the uncertainties by up to a factor of about 
1.5 in Figure 1 of that paper. Discussion in the earlier 
paper concentrated on uncertainties at Osft 2 = 0.019; 
at that baryon density, correction of the programming 
error reduces the estimated uncertainty on D/H by a 
factor of 1.6 and the estimated uncertainty on 7 Li/H by a 
factor of 1.2. Consequences of this error for cosmological 
implications are relatively unimportant, because the 
uncertainty estimates on our preliminary calculation are 
already smaller than those on astronomical observations. 

With regard to the nuclear data, we point out that 
any lingering problems with the nuclear database will 
probably have to be settled by new experiments. The 
low-energy p(n, j)d cross section is an important gap; 
the 3 He(c?,p) 4 He cross section is a less important gap, 
but it is clearly the case with the worst systematic 
problems in our nuclear database. Further reactions that 
contribute to lithium and deuterium yield uncertainties 
are listed in Sec. |y|. The importance of improving 
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various components of the nuclear database relative to 
improving the astronomical observations can be seen by 
examining Figs. [l^-[l7|. 

It is possible to improve on our work. We have explic- 
itly tried to work very directly with the nuclear inputs, 
and we have therefore avoided theoretical modelling, 
which is not essential in the face of such copious data. 
With more theoretical inputs — for example, radiative- 
capture models and multi-channel i?-matrix fits — one 
could include much more data, from higher energies and 
from scattering channels, and introduce some cross-talk 
between reaction channels. Solely in terms of gathering 
and processing data, this would be a much larger un- 
dertaking than what we have done. It could probably 
not be done piecemeal, because important theoretical 
parameters often have to be determined by examining 
data in several reaction and scattering channels, and over 
a wide range in energy. Depending on implementation, 
theoretical inputs may also present a more difficult 
problem in terms of propagating errors through to BBN 
reaction rates. 

Improvement may also be possible in our handling of 
correlated normalization errors. We have also avoided 
any fancier error estimation than applying quoted errors 
in a well-defined way. In particulary, we have not fit with 
floating normalizations because we were wary of altering 
the most important pieces of information for BBN. 
In place of floating-normalization methods, we have 
introduced a Monte Carlo method for making families of 
smooth curves to characterize data with normalizations 
that vary in the expected way. Although we do believe 
that this characterizes the uncertainties in the database 
in a reasonable way, it is worth noting that the individual 
curve fitted through each realization of the data does not 
recognize the presence of correlated errors. 

In conclusion, our results indicate that our method of 
coupling BBN error estimates to the nuclear data may be 
fruitful not only for providing a useful and unambiguous 
prescription for such error estimates, but also for making 
comparison between light-element abundances and BBN 
calculations more meaningful. While our prescription for 
handling the nuclear data is not unique, it is simple, 
repeatable, and direct. Such a method is needed for 
the program of "precision cosmology." We hope that 
our proposal will result at least in a more critical stance 
toward uncertainty estimates in BBN, and perhaps im- 
proved prescriptions that incorporate the better qualities 
of our method. 
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